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An improved e expansion in the d-dimensional (d > 2) stochastic theory of turbulence is con- 
structed at two-loop order which incorporates the effect of pole singularities at d — > 2 in coefficients 
of the e expansion of universal quantities. For a proper account of the effect of these singularities 
two different approaches to the renormalization of the powerlike correlation function of the random 
force are analyzed near two dimensions. By direct calculation it is shown that the approach based 
on the mere renormalization of the nonlocal correlation function leads to contradictions at two-loop 
order. On the other hand, a two-loop calculation in the renormalization scheme with the addition 
to the force correlation function of a local term to be renormalized instead of the nonlocal one yields 
consistent results in accordance with the UV renormalization theory. The latter renormalization 
prescription is used for the two-loop renormalization-group analysis amended with partial resum- 
mation of the pole singularities near two dimensions leading to a significant improvement of the 
agreement with experimental results for the Kolmogorov constant. 
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I. INTRODUCTION 



The renormalization-group (RG) method in the theory of turbulence allows to calculate various physical quantities 
— critical exponents and universal amplitude ratios — in the form of an expansion in a small parameter e. The real 
value of this parameter is not small, however, therefore justified doubts arise whether this method is of any use for 
^ \ acceptable numerical estimates of the quantities studied. Until recently, practical calculations were carried out only 
r**"* . in the simplest (one-loop) approximation and therefore it was not possible to assess how the next-to-leading terms 
of the expansion actually compare with the leading order at the real value of the parameter s. In Refs. [1] and [2] 
this problem was analyzed on the example of calculation of the skewness factor and the Kolmogorov constant in the 
inertial range. The calculation showed that the relative part of the two-loop correction is indeed large, of the order of 
100 % in the real space dimension d = 3. This contribution, however, rapidly decreases with the growth of d: already 
■ for d — 5 it yields only 30 % and in the limit d — > oo decreases to 10 %. On the contrary, when the space dimension 
decreases from d = 3tod = 2a drastic growth of the correction term is observed. 

Analysis of the dependence of the coefficients of the e expansion on the space dimension has revealed that this 
property is connected with the divergence of some graphs in the limit d — > 2, and the singularities in d — 2 = 2A 
accumulate with the order of the perturbation expansion. Contributions of these graphs turned out to give rise to the 
large value of the correction term also at d = 3. Thus, satisfactory quantitative results may be expected only after 
summing, at least approximately, the contributions of the most singular graphs at all orders of the e expansion. Such 
a summation has been carried out in Ref. [3] with the use of an additional renormalization and double expansion in e 
and A and with the result of a significant relative reduction of the correction term and improvement of the agreement 
with experiment. The calculation in Ref. [3] has been carried out in the two- loop approximation for both the usual 
e expansion and the double (e, A) expansion. These expansions were used as complementary to each other to arrive 
at the final result - an approach distinguishing Ref. [3] from Refs. [4, 5], in which the one-loop calculation in the (e, 
A) expansion was carried out. 

In Ref. [4], where the idea of the double expansion was first applied to the stochastic Navier-Stokes problem, the 
method used to eliminate the additional divergences was the natural at first sight multiplicative renormalization of the 
nonlocal correlation function of the random force (we shall further call this the nonlocal scheme). In Ref. [5] a different 
scheme of renormalization was adopted based on the general statement of the theory of UV renormalization that the 
counterterms are local. In the one-loop approximation, to which the authors of Refs. [4, 5] restricted themselves, it 
is possible to remove the divergences of the graphs both in the nonlocal renormalization scheme of Ref. [4] and in 
the local scheme of Ref. [5], so that at this level both approaches seem equally acceptable. This, however, is not so 
already in the next two-loop approximation. In this paper we show that the renormalization scheme of Ref. [4] is not 
consistent at two-loop order, because within it the two-loop contributions to the renormalization constants acquire 
an inacceptable dependence on the external wave vectors, whereas in the scheme of local renormalization [5] this does 
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not occur. The appearance of the latter scheme is more complicated and the reader unfamiliar with the delicacies of 
the RG approach might have doubts in the necessity to use it, all the more so because the argumentation in its favor 
requires a discussion of the so called A renormalization (similar to the shift of the critical temperature in the theory of 
critical phenomena). This paper is devoted to a detailed analysis of these issues. It also contains a technical account 
of the method which allowed to obtain the two- loop results announced in Rapid Communication [3]. 

This paper is organized as follows. In Sec. II wc remind basic features of the (field-theoretic) renormalization 
procedure and the subsequent asymptotic analysis in the two-loop approximation well above the problematic two 
dimensions. Sec. Ill is devoted to a detailed argument showing why the multiplicative nonlocal renormalization fails 
at the two-loop order of the double expansion. In Sec. IV the consistency of the local two-charge renormalization 
scheme is demonstrated by the results of the two-loop calculation in space dimensions d < 2 in which the technically 
simplest combined scheme of analytic and dimensional renormalization is unconditionally valid. The fairly technical 
issue of the possibility of analytic continuation of the results obtained for d < 2 to space dimensions above two is 
demonstrated in the Appendix. Renormalization-group equations are set up in Sec. V with the subsequent two-loop 
solution for asymptotic analysis in the inertial range. Details of the method of calculation of universal quantities in 
the improved e expansion are exposed in Sec. VI. Sec. VII contains discussion of the results and concluding remarks. 



II. RENORMALIZATION OF THE MODEL IN A FIXED SPACE DIMENSION d > 2 

The statistical model of the developed homogeneous isotropic turbulence of incompressible fluid is based on the 
stochastic Navier-Stokes equation 

Vm = v d 2 <Pi-diV + fi, v t = dt + (<pd). (i) 

Here, <fi(t, x) is the divergenceless velocity field, V(t, x) and fi(t, x) the pressure and the transverse random force per 
unit mass, respectively, v the kinematic viscosity. For the random force / a Gaussian distribution is assumed with 
zero mean and the correlation function 

</ i (t,x)/ j (t',x')) = A,(i,x;i',x') = j dkPii{k)d f (k) exp [ik(x-x')], (2) 

where Pij(k) = Sij — kikj/k 2 is the transverse projection operator, d the dimension of the coordinate space. For the 
function df(k) the following powcrlike form is adopted in the RG approach: 

d f (k)=D k 4 - d ~ 2E . (3) 

The quantity e > in Eq. (3) plays the role of a formal expansion parameter. The value corresponding to the physical 
model is e — 2, because for e — > 2, D ~ (2 — e) we arrive at dp(k) ~ <5(k) which corresponds to energy injection by 
infinitely large eddies. 

The stochastic problem (1) and (2) is equivalent to a quantum-field-theoretic model with a doubled set of transverse 
vector fields $ = {<p, <p'} and the action [6] 

= if/Diff/2 + <p'[-d t <p + v d 2 <p - (tpd)<p], (4) 

where D is the correlation function of the random force (2), and the necessary integrals over (i, x) and sums over 
vector indices are implied. Action (4) gives rise to the standard diagrammatic technique with the bare propagators 
whose (t, k) representation is of the form 

(<p(t)<p'(t'))o = 9(t - t') exp [-is k 2 (t - t')] , 

(<pV)o = 0, (5) 

(w)o = |^exp[-i,ofc 2 |*-*'|] , 

where the common factor -Pjj(k) has been omitted for simplicity. The interaction in Eq. (4) brings about the three- 
point vertex —ip'(ipd)ip = ip^Vijsipjips/2 with the vertex factor Vij S = i(kjSi s + k s dij), where k is the wave vector of 
the field tp'. The expansion parameter of the perturbation theory is the coupling constant go = Dq/vq. 

Model (4) is logarithmic (i.e. the coupling constant g is dimensionless) at e = 0. In the analytic renormalization 
scheme adopted here the UV divergences have the form of the poles in e in the correlation functions of the field $ = 
{ip, ip'}. Dimensional analysis (power counting) shows that for d > 2 superficial UV divergences can be present only 
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in the one-particle-irreducible (1PI) functions r v > v and Y v i vv . These divergences may be removed by counterterms 
of the form 

<p'd 2 <p, <p'dt<p, <p'(<pd)<p (6) 

in the action. Due to symmetry reasons, however, in model (4) only one counterterm of all allowed by the dimensional 
analysis is actually generated. First, the spatial derivative acting on the field ip in the interaction term of action (4) 
can be transferred to the field <p' with the use of integration by parts. This means that the counterterms to the 1PI 
functions must contain at least one spatial derivative, so that the structure (p'd t (p cannot possibly appear. Second, 
from the Galilean symmetry of action (4) it follows that the last two structures of Eq. (6) can be brought about as 
counterterms only in the invariant combination ip'Vt<p with the Lagrangian derivative Vt = <9t + (fd) from Eq. (1). 
This excludes also the structure ip'(ipd)(p. Thus, in the generic case we are left with a single counterterm of the form 
ip'd 2 ip. In the special case d = 2, however, a new UV divergence appears in the 1PI function Y v i v >. 

Consider the renormalization of model (4) in the two- loop approximation in d > 2. In this case the only counterterm 
required is ip'd 2 ip which is generated by multiplicative renormalization of the viscosity in the corresponding term 
of action (4). We shall use the scheme of minimal subtractions (MS) in which the renormalization constants are 
determined by the relations 

,3 _ „,,2e,,3 



v Q = vZ v , D = g Vo = g\i V 

7j gi 



go = g^Z g , Z g = Z-\ (7) 



Here, \i is the scale-setting parameter (the reference mass) in the MS scheme, v is the renormalized viscosity and g the 
dimensionless renormalized charge. The only independent renormalization constant in Eq. (7) is that of the viscosity 
Z v . The amplitude of the correlation function of the random force D a is not renormalized, because no counterterm 
of the form (p'(p' in action (4) is necessary. This leads to the relation between the renormalization constants of the 
charge and viscosity indicated in Eq. (7). 

In the MS scheme the renormalization constants are constructed as Laurent series in e of the form "1 + 
J2n>i a " £ ~ n " ■ I n particular, 
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Z v = l + u^f+u 2 (^ + ^i) + ... = l + ^ u ™^a„ fc ^ fc , (8) 

where 



71=1 fe=l 



= gSd , _ _Sd_ „ = 2^/ 2 
U - 32 ' d ~ (2^) d ' d ~ T{d/2) ' 1 ' 

and the coefficients a n k depend only on d. Here Sd is the surface area of the unit sphere in d-dimensional space and 
r is Eulcr's Gamma function. 

We shall determine the constant Z v from the requirement that the 1PI correlation function T^,^ at zero frequency 
(ijj = 0) is UV finite, i.e. finite at e — > when expressed as a function of the renormalized variables v and g determined 
by relations (7). With respect to vector indices the function T v i v is proportional to the transverse projector Py (p), 
where p is the external wave vector. In the following we shall deal with the scalar coefficient of this projector 
obtained by the contraction of the indices i and j and division by TrP = d — 1. In terms of the bare parameters v$ 
and Do = ffo^o this scalar coefficient at lo — assumes the form — vop 2 + sum of contributions of the n-loop graphs, 
each of which contains n pieces of (ff)o lines (5) and, correspondingly, the factor . Thus, in view of dimensional 
arguments 



Tr r^'m \oj=o 2 
d-1 = "° P 



E( DpSd \" (n) 



(10) 



with dimensionless coefficients 7^ which only depend on d and e. The factors 32 and Sd in Eq. (10) have been 
introduced for convenience. To obtain the renormalized function T V ' V the parameters Do and vo in Eq. (10) have 
to be expressed in terms of u, g and fi according to definitions (7), which leaves the coefficients 7^ intact. It is 
convenient to divide the result by vp 2 to arrive at the dimensionless quantity 

^^iy = + - 2e 1% ^ 2 + (- 2£ ) 2 ^ ^ + - (") 



4 



with u from Eq. (9) and s = fi/p. 

The renormalization constant Z v is determined from the condition of cancellation of the poles in e in relation (11). 

In the coefficient j^) there is a simple pole ~ 1/e, whereas j^) contains poles ~ 1/s and ~ 1/e 2 etc. For the 
two-loop calculation of Z v the following contributions are needed: 



y (2) - 



+ 5 + 



13 



+ ... 



(12) 
(13) 



where the ellipsis stands for irrelevant corrections 0(e) in 7^ and 0(1) in 7^. 

Denoting the contribution of the order u n ~ <?" to the renormalization constant (8) by zi n \ from the condition of 
cancellation of the divergences (poles in e) in Eq. (11) we infer 



zP=c e 



(14) 
(15) 



where C e stands for the operation of extraction of the UV-divergent part, which here consists of poles in e. 

When relation (12) is substituted in Eq. (14) the UV-finite term B does not contribute and the coefficient s 2e 
1 + 2elogs + ... may be replaced by the unity. As a result we obtain 



uA 



(16) 



Substituting this expression together with relations (12) and (13) in Eq. (15) we find 



zip 



u 2 s 4s 



-z + — - 2u 2 s 2e - — 



B 



(17) 



In the terms <~ 1/e we may replace s ne — » 1, whereas in contributions ~ 1/e 2 also the second term in the expansion 
g ne _ 2 7 nelogs 7 ... must be retained which gives rise to a contribution of the form e" 1 logs = e _1 log(^/p) in 
(2) 

zr. The presence of such a term in Z v is inacceptable, because renormalization constants must not contain any 
wave-number dependence by their very definition. The condition of vanishing of the term <~ e _1 logs in (17) is 

C = A 2 



for the coefficients of relations (12) and (13). 

The recent two-loop calculation [1] confirms that relation (18) holds. Substituting it in Eq. (17) we obtain 



ZiP 



^2 



D - 2AB 



(18) 



(19) 



The one-loop coefficient A in Eqs. (16), (17) and (19) has been known for quite a while: 

A __Md-i) 

A ~ d + 2 ■ 

For the nontrivial next-to-leading coefficients D and B in Eq. (19) integral representations readily calculable for any 
given d have been obtained in Ref. [1]. 

That condition (18) holds thus imposing on Z v cancellation of the contributions ~ logs is not a coincidence, 
but a consequence of general principles of the theory of UV renormalization. The most important of them is the 
requirement that all counterterms must be local in space (i.e. polynomial in wave vectors). In model (4) this is so, 
because the counterterm giving rise to the renormalization of the parameter vq has the form of vp 2 multiplied by a 
wave-number-independent coefficient, i.e. a polynomial function in p. Therefore in this model all consequences of the 
general conjectures of the theory of UV renormalization must hold, in particular, independence of the renormalization 
constants of wave numbers to all orders in the perturbation theory as well as the critical scaling due to the RG 
equations with the e-dependent critical dimensions of the velocity field ip and the frequency u) (more details in Sec.V): 



A ip = l- 2e/3, A w = 2 - 2e/3. 



(20) 
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These are exact relations without any corrections of higher order in e. They are a consequence of connection (7) 
between the renormalization constants Z g and Z v which, in turn, follows from the absence of renormalization of the 
nonlocal contribution with the correlation function of the random force in action (4) . At the real value e — 2 quantities 
(20) assume the Kolmogorov values 

A v = -1/3, A w = 2/3. (21) 

Condition (18) ensuring independence of the renormalization constant of the wave number in the MS scheme may 
appear in a different form in other renormalization schemes. We will illustrate this point on the example of the 
scheme with the "normalization point" (NP). In this approach the renormalization constant Z v is calculated from the 
normalization condition for the 1PI Green function 



TrT 



<p'<p \iO — 



vp 2 {d-l) 



= -1 (22) 



in contrast to the cancellation of poles in e in expression (11) in the MS scheme. Then instead of Eqs. (14) and (15) 
we obtain 



4 1} = «7& , 



~ 2(7^) a , (23) 



and after substitution of Z v from Eq. (23) expression (11) assumes the form 



vp 2 (d-l) "f'f 



= -1 + «7#l (* 2e - 1) + u 2 j { X( S is - 1) - 2( 7 « )V £ - 1) + 0(u 3 ). (24) 



In the NP scheme the renormalization constant Z v does not depend on s = fi/p due to the very definition, but 
cancellation of poles in e is not obvious in Eq. (11). In the two-loop approximation (24) with account of expressions 
(12) and (13) these poles appear in the form ~ it 2 £ -1 logs in several contributions, and the condition of their mutual 
cancellation is the same relation (18) which ensured the cancellation of the "bad" contributions <~ u 2 e _1 logs in Z v 
in the MS scheme. As it was previously explained, fulfilment of condition (18) is guaranteed by general theorems of 
the theory of UV renormalization with local counterterms. 

The MS and NP schemes differ by a finite renormalization of the parameters g and v, therefore all objective physical 
quantities, in particular, critical dimensions (20), calculated in these schemes coincide. 

Critical dimensions (20) do not depend on d and thus for them the problem of singularities in the limit d —> 2 
mentioned in Sec. I is not relevant. There are, however, other important physical quantities such as the skewness 
factor, Kolmogorov constant, critical dimensions of various composite operators to which this problem persists. It 
is important that for these quantities the problem of anomalous scaling is absent, which cannot be treated in the 
framework of the model with massless injection (3) lacking a dimensional parameter to account for the external scale 
of turbulence. 

For such quantities, contrary to Eq. (20), the solutions contain full series of the form 

oo 

R{e,d) = Y,Md)e k , (25) 

fe=0 

and the coefficients Rk(d) in the limit d — > 2 reveal singular behavior of the type ~ (d — 2)~ k ~ A~ k (2 A = d — 2) 
giving rise to the growth of the relative part of the correction terms at d — > 2. The effect of these is fairly discernable 
also at the real value d — 3, hence the natural desire to sum up contributions of the form (e/A) k at all orders of the e 
expansion (25). This may be done with the aid of the double (e, A) expansion [4, 5]. The idea of such an "improved 
e expansion" with the use of the local renormalization scheme [5] was explained in our Rapid Communication [3], 
where many important subtleties and details of calculations were, however, not reflected due to lack of space. In the 
present paper we give a detailed exposition and start from the proof of inconsistency of the renormalization scheme 
proposed in Ref. [4]. 



III. CONSTRUCTION OF THE DOUBLE (e, A) EXPANSION. PROOF OF THE INCONSISTENCY OF 
THE NONLOCAL RENORMALIZATION [4] IN THE TWO-LOOP APPROXIMATION 

Model (4) is logarithmic (i.e. the bare coupling constant go is dimensionless) at e — in function (3) in arbitrary 
space dimension d. In a fixed dimension d > 2 the value e = 2 corresponds to the "real problem". Calculations 
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in the framework of the e expansion have a rigorous meaning only in the vicinity of e = 0, whereas continuation of 
the results to the "real" value e — 2 is always understood as an extrapolation. In the scheme applicable for d > 2 
reviewed in Sec. II this extrapolation corresponds to the continuation along the vertical ray from the point (d, e = 0) 
to the point (d, e — 2) in the (d, e) plane. The same final point may be reached along a ray from any starting point 
(do d, e — 0) at which the model is logarithmic as well. The extrapolation along the ray starting from the origin 
(do = 2, e = 0) is, however, singled out, because at d = 2 in model (4) an additional UV divergence (absent at d > 2) 
occurs in the 1PI function T v i v i. On such a ray we put 

rf = 2 + 2A, A/e = ( = const. (26) 

The parameters e and A are considered small of the same order and their ratio A/e = (a fixed constant [( = 1/4 in 
the extrapolation to the point (d = 3, e = 2)]. 

Extraction of contributions of the order e m with A/e — const corresponds to the account of all contributions of the 
form e m (e/A) n with any n = 0, 1, 2... and m + n = k in Eq. (25). Thus the use of the (e, A) expansion in such a 
form is directly related to the problem of the account of the singularities at A — > pointed out in the discussion of 
relation (25). 

It is worth emphasizing that the very process of extrapolation along a ray from the starting point (d = 2, e = 0) 
is inapplicable to description of two-dimensional turbulence in which the physics is totally different from the three- 
dimensional problem due to the appearance of the inverse energy cascade [7] . In Fig. 1 we have plotted the borderline 
curve B AC between the direct (normal) and inverse energy cascades obtained in Ref. [8] . The starting point of the 
extrapolation for the two-dimensional case (d = 2, s = 0) lies in the region of the direct cascade, whereas the final 
point (d = 2, e = 2) in the region of the inverse cascade. Thus the ray connecting these points intersects the borderline 

the curve BAC - so that the extrapolation becomes impossible. However, the ray connecting the starting point 
(d = 2, e = 0) and a final point like (d = 3, e = 2) lies completely in the region of the direct cascade, therefore on 
such a ray the problem of the change of the cascade pattern does not arise. The rightmost point of the region of the 
inverse cascade (point A on Fig. 1) has the coordinate dA — 2.06 [8]. In the preceding discussion of the extrapolation 
along the vertical ray from the point (d, e = 0) to the point (d, e = 2) at d > 2, it should have been noted that the 
condition is not simply d > 2, but d > d^ — 2.06. From the practical point of view this is irrelevant, because we are 
interested in the space dimension d = 3. 

The idea of the double (e, A) expansion together with the extrapolation along the ray A ~ e of relation (26) in the 
context of the present problem was first put forward in Ref. [4] . The UV divergences are present not only in the 1PI 
function T v i v but also in and appear in the form of poles in the parameters e and A and linear combinations 

thereof, or, cquivalently, as poles in e with the fixed ratio A/e = ( = const. To remove the additional divergences 
from the graphs of the 1PI function T v i v i renormalization of the amplitude Dq in the nonlocal correlation function of 
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the random force (2) and (3) was used in Ref. [4], i.e. relations (7) between bare and renormalized parameters were 
replaced by 



vq = vZ v , D = ,g ^o = 5M e v A Z D 



g Q = gn 2 "Z g , Z a Z% = Z D (27) 



with a new renormalization constant Zd which does not have an analog in Eq. (7). 

It should be noted that the introduction of the additional constant Zu breaks the last connection in Eq. (7) and 
its consequences (20). Therefore, the author of Ref. [4] has put forward the conjecture that in the scheme of the 
double (e, A) expansion at the real value e = 2 the velocity field ip and the frequency uu have dimensions with values 
different from the Kolmogorov values (21). This is, of course, true, if renormalization relations (27) are used. We shall 
further show, however, that the renormalization scheme of Ref. [4] with relations (27) is not internally consistent. 
This is not obvious in the one-loop approximation, to which the author of Ref. [4] restricted himself, but becomes 
apparent already in the next two-loop approximation. In Ref. [5] another scheme of construction of the double (e, 
A) expansion was put forward in which the last equality in Eq. (7) together with its consequences (20) and (21) are 
preserved. We shall deal with this approach in Sec. IV. 

The main goal of this section is to prove that the scheme of multiplicative renormalization (27) contains intrinsic 
contradictions. To this end, consider representations similar to (11) for the 1PI functions T^^ and T^^. According 
to Eq. (27) the amplitude D in Eq. (10) now acquires the additional factor Zjj, therefore instead of relation (11) 
we now obtain 

= + ^ Z " Zd + (« s2£ ) 2 Z ^ Z n + - ' ( 28 ) 
The analogous relation for the 1PI function r v ' v ' is 

g^^epi-d^e^ — — Z D + US 7^, Z v Z D + (us ) 7^, Z v Z D + ... . (29) 

The expansion parameter is u = gSd from Eq. (9). In Sec. II the quantity d was considered a fixed parameter and 
therefore it was possible to treat Sd as a simple normalization factor. Here, d is determined by the relation (26) 
and in calculations within the usual MS scheme the quantity Sd should be expanded in the small parameter A ~ e. 
Following Ref. [4], we shall use the modified scheme MS (see, e.g., Ref. [9]), in which the quantity Sd is treated as a 
whole and not expanded in A. It is well known that the choice of scheme is not reflected in any physically significant 
results. 

The constants Z are sought as series of form (8) and determined from the condition of cancellation of the UV 
divergences (poles in e with A/e = const) in relations (28) and (29). Denoting by Z^ n > the contribution of order 
u n ~ g n in any of these constants we arrive at expressions similar to (14) and (15): at the first order in u ~ g 

Z$=-C e {us^%}, (30) 

and at the second order 



4 2) = £s {u's^X + us Hl i Z n 2^>] } , (31) 

4 2) = {-u^X' + us Hl' ^ - ^ ■ (32) 
For calculation in the two-loop approximation the following contributions are needed 

+ B + 7&-y + * + ". (33) 

1&-? + 7 + ". 7&, = f + 7+<"- (34) 

These are analogs of relations (12) and (13) with different coefficients, however, which now may depend on the ratio 
A/e = 
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Substituting expressions (33) in Eq. (30) we find the one-loop contributions to the renormalization constants: 

^ = ^, Z$ = -^. (35) 
One-loop calculation yields the following values (first obtained in Rcf. [4]): 

A = -l, A' = ^. (36) 

In the one- loop approximation there are no problems with logs in the constants Z, so that the multiplicative renor- 
malization (27) appears quite acceptable. 

Consider now two-loop contributions (31) and (32). Taking into account the already known one-loop expressions 
(35) we obtain 

Z» = C, LV + us- (- + b) (-— - — ) 1 , (37, 



e l \ e e 



^M-»v<(^)w.(f + *.)(^ + ^)}. 

The condition of cancellation of the contributions ~ e _1 logs in Eq. (37) is 

AC + 2A{-A' -2A) = 0, (39) 

and in Eq. (38) analogously 

-4C' + 2A'(3A + 2A r ) =0. (40) 
Our two-loop calculation of the coefficients C and C yields 

1 



G = \- 



2(2 + 0' 



C = — . (41) 

(2 + C)(3 + C) (3 + 1 ' 

Substitution in relations (39) and (40) of the calculated quantities (36) and (41) readily shows that condition (39) 
is satisfied, whereas (40) is not. This means that in there is no "bad" contribution <~ e 1 logs = e 1 log(/z/p), 

(2) 

while in Z D there is such a term: 

2(1 + 0(4 + 30 , 

(2 + 2 (3 + " gWP) ' ( ) 

whose coefficient is the expression of the left-hand-side of Eq. (40). 

Thus within the renormalization scheme of Ref. [4] according to relations (27) a dependence on the external wave 
numbers through logs = log(/x/p) appears in the renormalization constants, which is completely inacceptable by the 
very definition of the renormalization constants. It is not difficult to understand the reason of this: in scheme (27) 
there is a violation of a fundamental principle of the general theory of UV renormalization - requirement that all 
counterterms must be local (polynomial functions of external wave vectors) [9] . The introduction of the coefficient 
Zd at the term ~ ip' 'Dip' in action (4) with the nonlocal injection function (3) is tantamount to introduction of a 
nonlocal counterterm with the structure p A - d - 2e . This feature takes the scheme discussed beyond the framework of 
the standard theory of UV renormalization with such unpleasant consequences as the appearance of the (inacceptable) 
dependence on wave numbers in the renormalization constants. This general line of argument motivated the authors 
of Rcf. [5] to change the scheme of (e, A) renormalization to conform to the requirement of the polynomial in wave 
vectors form of all the counterterms (localness) , although in the one-loop calculation of Ref. [4] the inconsistency of 
the scheme proposed there does not show explicitly. 

( 2) 

It might be suggested to change relation (32) to exclude the wave-number-dependent contribution (42) from Z K D ' . 
Eq. (32) was obtained, however, from the requirement that in the two-loop approximation all UV divergences - poles 
in £ - were removed from the renormalized 1PI function T v > v /, so that any change of the form of from (32) would 
lead to the appearance of poles in e in the renormalized function T^/^i. 
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The persistent opponent might say:" Who cares, I am not interested in the two- loop approximation, I am completely 
happy with the one-loop accuracy, where there are no problems." Here, the objection would be that elimination of 
UV divergences (poles in e) to all orders in perturbation theory is not a caprice but a compelling necessity. If such 
poles are left, then there is no guarantee that results obtained at the lowest order of perturbation theory do not 
acquire corrections of the same order from the higher-order terms not accounted for (in fact, there is conviction in 
the opposite), i.e. lowest-order calculations become completely unreliable. Therefore, in particular, the conclusion of 
Ref. [4] that relations (20) are violated in the (e, A) scheme is not correct; in the consistent renormalization scheme 
these relations continue to hold [5]. 

Let us briefly discuss the possibility to carry out a nonlocal renormalization in the NP scheme. Relations (27) are 
preserved in this case, whereas the two independent renormalization constants Z v and Zd 2 are to be determined from 
the following normalization conditions at p = [i for the 1PI functions (28) and T^t^t (29): 

Tr r^'^l^o 
vp 2 [d- 1) 

5 z/V e p 4 - d - 2e (d- 1) 

The problem of dependence of the renormalization constants on wave number is absent in such a setup. However, it 
may be readily checked that conditions (39) and (40) remain necessary to ensure absence of UV-divergent (at e — > 0) 
contributions <~ u 2 e _1 log(yu/p) in the renormalized Green functions r v < v and T v > v i for arbitrary values of the wave 
number p. 

(2) 

In conclusion, let us point out that the "bad" contribution (42) in Z D ' vanishes at ( = —1, i.e. at A = — e in 
Eq. (26). Then d = 2 + 2A = 2 — 2e and energy injection (3) becomes local: df <~ p 4 - d ~ 2e — p 2 (such a model 
was considered in Ref. [10]). In this case the multiplicative renormalization (27) conforms to the requirement of 
local counterterms and the corresponding constants Z do not contain any dependence on log s in accordance with the 
general theory. 

IV. CONSTRUCTION OF THE (e, A) EXPANSION IN THE TWO-CHARGE MODEL WITH LOCAL 
COUNTERTERMS. TWO-LOOP CALCULATION OF THE RENORMALIZATION CONSTANTS 

In the preceding section it was shown that in the (s ,A) scheme (26) the multiplicative renormalization [4] of 
the amplitude D in Eq. (3) is not acceptable. The reason is that the counterterm with structure (3) is nonlocal 

_ fc 4-d-2s = fc 2-2A-2e Qn rays (26). 

Guided by the general theory of the UV renormalization, the authors of Ref. [5] put forward another scheme, in 
which a local counterterm <~ k 2 instead of the nonlocal one ~ /j2-2A-2e j g usec j ^ aDSOr b singularities from the graphs 
of the 1PI function T v ' v '. This corresponds to addition of the term ~ (p'd 2 ip' to the action functional. In functional 
(4) with the correlation function D from Eqs. (2) and (3) there is no such term, so that upon the addition of the 
term ~ ip'd 2 ip' the renormalization ceases to be multiplicative. This would be unessential, if our only goal was the 
elimination of divergences from Green's functions which is quite possible by a non-multiplicative renormalization. For 
the use of the standard technique of the RG multiplicative renormalization is, however, necessary. This is why the 
authors of Ref. [5] proposed to consider a two-charge model in which to function (3) ~ fc 4 ~ d ~ 2e = fc 2 -2A-2e ^ ^ erm 
~ k 2 is added at the outset with an independent coefficient: 

df {k) = D 1Q k 2 - 2A - 2 ^ + D 2Q k 2 = g W ^ fc 2-2A-2s + g2QU 3 k 2 (44) 

Here, the amplitude D n of Eq. (3) is denoted by D w . The parameters gio and 520 introduced in Eq. (44) play the 
role of two independent bare charges. 

The contribution with D20 in relation (44) corresponds to thermal fluctuations. A model with this term only has 
been analyzed earlier in Ref. [10]. In the theory of turbulence D 2 o = should be considered the "real value" of this 
parameter, since only the first term in Eq. (44) at e = 2 reproduces the realistic for the theory of turbulence pumping 
of energy by large-scale eddies. It will be shown below that vanishing of the bare parameter g 2 o = £>20^cT 3 = does 
not imply vanishing of the corresponding renormalized parameter g 2 , so that in terms of renormalized parameters 
function (44) gives rise to a two-charge model. 

The unrenormalized action is, as before, functional (4), but now with the injection function (44) instead of (3) in 
correlation function (2). In the adopted shorthand notation 

= ^(D w k 2 - 2A - 2s + £> 20 fcV + <p'[-9t<P + vod 2 tp - (tpd)<p] . (45) 




= 1. 



(43) 
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The propagators {<p<p')o and (tp'tp')o corresponding to action (45) maintain earlier form (5), whereas (<p<p)o is replaced 

by 



(D 



j,2-2A-2e 
10 K 



D 20 k 2 ) 



2v k 2 



exp [— u k 2 \t — t'\ 



(46) 



We are interested in the region e > and A > in Eq. (26). In this region in model (45) the additional problem of 
"A divergences" arises which was absent in model (4) with injection function (3). Let us explain this in more detail. 
Wave-vector integrals - with the shorthand notation J dk . . . - corresponding to the 1PI graphs discussed always 
reduce to "nearly logarithmic" ones in the present set of models. Their deviation from logarithmicity shows in the 
form of factors of the type k a with a small exponent a — 2mA — 2ns, where n and m are nonnegative integers. The 
exponent a is the wave-number dimension of the wave-vector integrals obtained upon all time integrations and may 
be calculated by the following simple rule: each loop integral over wave vectors contributes a term 2A to a, the term 
with D w in Eq. (46) yields the contribution — 2e — 2A, but the term with D 20 does not affect a at all. Thus it may 
readily be seen that if only the nonlocal term with Diq is left in Eq. (44) [i.e. if we return to model (3)], then all the 
exponents a in graphs of T v i v and Y^i^i at s > and A > become negative. All the integrals in the limit k — > oo 
converge, they may be carried out over the whole wave-vector space and the divergences show as poles in e, A and 
their linear combinations. 

However, in the model with injection (44) - due to the presence of the second term with D 2 q - at A > wave- vector 
integrals appear with a > 0. They diverge in the limit k — > 00 and thus require an UV cutoff A. As examples we 
quote the values of a in the graphs of interest for us. In the one-loop graphs of T v ' v : a — —2s, 2 A; in the two-loop 
graphs: a = — 4e, — 2e + 2 A, 4A; in the one-loop graphs of r v ' v i: a = — As — 2 A, — 2s, 2A; in the two-loop graphs: 
a = -6s - 2A, -As, -2s + 2A, 4A. 

Thus, in the two-charge model (44) at A > (e > is always implied) some integrals have the A divergence at 
large k. To remove these divergences an additional procedure of A rcnormalization procedure is needed which we shall 
discuss in the Appendix. At the moment the important point is that after the A renormalization the limit A — ► 00 
may be taken with the result that divergences appear only in the form of poles in s, A and their linear combinations. 
The same poles may be found within the " formal scheme" , where all integrals are understood as analytic continuation 
on the parameters s and A from the region, where there are no A divergences. 

In our case this is the region of s > and small (compared with e) negative A < (i.e. d < 2). In this section 
we shall consider results obtained in the framework of this " formal scheme" . There is no UV-cutoff parameter A in 
this scheme, but the divergences appear in the form of poles in s with A/s = const. The goal of the renormalization 
is removal of these poles. In the Appendix it will be shown that the results obtained this way coincide with those 
obtained at A > after the A renormalization and subsequent limit A — > 00. 

The relations of multiplicative renormalization in the formal scheme are 



9w"o 



D 20 
vq = vZ v , 



gin 25 " 3 , 



320^0 



92V 



-2A,,3 



v A Z 



D 2 , 



9w 
920 



9iV 2£ Z gi , 
g 2 u- 2A Z g2 , 



(47) 



with two independent renormalization constants for the viscosity vq and the amplitude -D20; the amplitude Diq of the 
nonlocal correlation function of the random force is not renormalized. The renormalization constants Z v and Zu 2 arc 
found from the condition that the 1PI functions | w=0 an d Y v ' v ' \ u —q are UV finite (i.e. with A/s = const there 
are no poles in e). The dimensionless expansion parameters of the perturbation theory for these quantities are 



DwSd 



32*gp 2e ' 

with Sd from Eq. (9). Instead of relation (10) we now have 



a 2 



D 2 oSd 



32i/gp 



-2A 



(48) 



TrT 



<£>'<£> \UJ — 9 

in— = U0P 



-1+ E «i la 2 2 7i"r 2) 



l>0,»2>0. 
tl + n 2 >l 



(49) 



and the analogous expression for T 



Tr r^/^y |w=o 
1^1 



= D w p 



2-2A-2e 



D 2 qp 



1+ ]T ar«2 2 7?r 2) 



(50) 
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In terms of the renormalized variables relations (49) and (50) yield for the reduced dimcnsionless functions the 
following representations 

TrTyV L=0 _ 7 > 7 ST «l « 2 («i,n 2 ) /ri\ 

F V ' ji 1 >0,n 2 >0, 

„ 1 +„ 2 >1 



Tr r^/^j/ o 
(d — l)g 2 v 3 u~ 2A p 2 u 2 



(52) 



■>l>0,n 2 >-l, 
»1+«2>1 



where the expansion parameters a\ and «2 from Eq. (48) are expressed through the renormalized parameters according 
to relations (47): 



-2A r 



(53) 



Here, u\ = giSd/32, u 2 — g 2 Sd/32 and s = fi/p. Dependence on e of the coefficient functions 7^'" 2; and 7^'/' 
in Eqs. (51) and (52) is determined by relations of the form of Eqs. (33) and (34), in which the £ = A/e-depcndcnt 
coefficients A, B, C, D, A', B' , C and D' now acquire subscripts corresponding to the superscripts (711,712) of the 
quantities 'yC™ 1 -™ 2 ) . In the one-loop approximation the following analogs of relations (33) are needed 



(»,fe) _ 



A 



i,k 



+ B i k , 



A' 



i,k 



+ B 'i,k . 



with the index sets (i, k) = {(1, 0), (0, 1)} for 7^ and the sets (i, k) = {(2,-1), (1, 0), (0, 1)} for 
approximation the following analogs of relations (34) have to be included 



(54) 

In the two-loop 



(i,k) _ Cj.fc Dj,k 
e 



"V<P £ 2 



(*>fe) 
7 vV 



C ik 



+ 



(55) 

with the sets (i,k) = {(2, 0), (1, 1), (0, 2)} for 7^ and (i,k) = {(3, -1), (2, 0), (1, 1), (0, 2)} for 7 vV . For calculation 
of the constants Z in the one- loop approximation, which was carried out in Ref. [5], only constants A and A' from 
Eq. (54) are needed. For our two-loop calculation all constants in Eqs. (54) and (55) are necessary. 

The constants Z v and Zd 2 are determined from the condition of cancellation of all UV divergences ( poles in e with 
A/e = const) in Eqs. (51) and (52). Denoting by Z^ the contribution of order u n ~ g n with respect to the set of 
charges Ui and u 2 to any constant we obtain at the first order in u ~ g 



7(1) 

J D 2 



-2A (0,1) 
U 2S 7^ V 



U l 4e+2A (2,-1) , 2e (1,0) , -2A (0,1) 



(56) 



and at the second order 



2 4e (2,0) . 
"I s Xp' v ' +U 1 U 2 S 



2 s 



-2A (1,1) 
'<P'<P 



+ UnS 



-4A (0,2) 



-2Z« 



+ u 2 s 



-2A (0,1) 



V -D 2 



2Z. 



?>]}, (57) 



7(2) 

J r> 2 



l U2 



6£+2A (3,-1) 
* Itp'tp' 



2„4e (2,0) 



+ !fl s 4e+2A (2,-1) 
«2 



/ (ft (ft 



MP 



2e-2A„,(1,1) 



2„-4A (0,2) 



tp'tpi 



U 9 S 



Z% 3Z« 



-2A (0,1) 

+ u 2s i y v ;J 



2Zg - 3Z« 



• (58) 



Substituting expressions (54) and (56) we find the one-loop contributions to the renormalization constants 
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= - ( Ul A h0 + u 2 A 0A ) , Z%1 = -- f — ^2,-1 + «i^i,o + u^.i) 
The coefficients A and A' here have been calculated in Ref. [5]. In our notation 

^i,o = — 1) = t , ^2,-1 



(59) 



2 + C 



4' — 9 
^1,0 — z ' 



A' — - 

-^0,1 — 



c 



(60) 



In the present work we have carried out the two-loop calculation and determined the coefficients B and B' in Eq. 
(54) together with C, C , D and D' in Eq. (55). Let us quote the coefficients C and C necessary at the moment: 



C2,o — 1 

^3,-1 — 



1 



r 2 r 1 



c 2 , Q = -i + 



2(2 + 

2 3,1 

(2 + C)(3 + C) ~ 3 + C ' C '°- 2 = _ 2?J ' 

1 3 4 j 



(61) 



2 + C 2( 



C(i-C) i-C 



To check the cancellation of the "bad" terms <~ e _1 logs in Eqs. (57) and (58) only terms ~ 1/e 2 from them arc 
needed. They are determined in Eq. (55) by the coefficients C and C from Eq. (61) and in the contributions with A 
and A' from Eqs. (54), (59) and (60). Substitution shows that all contributions with e _1 logs in Eqs. (57) and (58) 
cancel as required. 

A specific feature of the renormalization constant Zd 2 is that it contains terms <~ l/ti 2 [see Eq. (59)]. When 
such a Zd 2 is substituted in renormalization relations (47) in the expression for D 20 terms independent of u 2 appear 
(generation terms): 



U2 - - {u\ A' 2 _j + utu 2 A' 10 + U^O.l) + ■ 



(62) 



Due to such terms the condition D 2 q — does not lead to the trivial conclusion u 2 = 0, i.e. a nonvanishing value of 
the renormalizcd charge corresponds even to the zero (real) value of the bare charge. 

The ellipsis in Eq. (62) stands for contributions of the two-loop order and higher, which contain terms ~ 
with n > 3. In the region u <~ e (where the fixed point u* of the RG lies, see Sec. V) they are of the same order 
in e as the explicitly quoted one-loop contribution in Eq. (62). Therefore, to determine the connection between the 
charges u\ and u 2 imposed by the condition D20 = (i.e. Zd 2 = 0) the two-loop calculation of the constants Z is not 
sufficient. This is unimportant, however, in the following, because in the RG analysis of Sec. V the charges u\ and 
u 2 are considered independent parameters. 

We shall not quote here the fairly cumbersome expressions obtained by us for the constants B, B' , D and D' in Eqs. 
(54) and (55). Instead, we quote the two-loop expressions for the renormalization constants Z v and Zd 2 obtained 
with the use of them and relations (57) - (61) in the MS scheme (a detailed account of the method of calculation can 
be found in Ref. [1]): 



Z v = 1 



Ml U2 



4C + 3 2C+1 



(2 + 0- 



(e 2 



Ml 



5C + 3 



+ 



e(l-C) (l-C)e 



U\U 2 



I— 1 


1 


1 


2 1 


2 




C 2 e 2 _ 


M2 + - 

e 



Ul 



U\U 2 

: (i-0 



u 2 

c 



R, (63) 



Ul 



J D 2 



e«2(2 + C) 
34C+ 19 + 6C 2 

(2 + 0* 



2ui u 2 



C (13 + 19Q 
2(3 + (2 + 0e 



+ 



(C + 4) (2C+1) 



C(2 + * 2 



V- 



13 



2C + 1 



(3 + (2 + 0^ 2 
b 31 C 2(4( + 



U\ 3 U 2 ~ 



1) 



e(l-C) (l-0Ce 2 



UlU 2 



+ 



Ui 



J w 2 (3+ 



3 ui +6 



U\U 2 

(W) 



3 

ul 

c 



1 



C 2 e 2 
(R- 



u 2 



1), (64) 
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where 



R 



-0.168. 



This number has been obtained by a computer calculation of a relatively simple but cumbersome twofold integral, 
through which all the nontrivial two-loop contributions in Eqs. (63) and (64) are expressed. 

In conclusion, we quote the analogs of relations (63) and (64) in the NP scheme. In this scheme the constants Z 
are determined instead of Eq. (43) by normalization conditions 



Tr r^'cp \tj— o 



vp 2 (d- 1) 



-1, 



gv 3 H 2s p 4 - d - 2 c(d~l) 



Ui 

u 2 



+ 1, 



(65) 



with T v i v from Eq. (51) and T^i^i from Eq. (52). From here the renormalization constants follow in the form: 



Z v = l + 



~-§C + c) ": ' (77 i! 1 ' 



+ 
+ 



(C-i)e 2 



2c- 3C 



6 + 2c- - 



R- 
R 



c-i 



uiu 2 



Ul 



2( 2 e 2 



R 



U2 



(66) 



Zd 2 =1 



(2 + 0* 



7 5 

2 + (2 + C) 



U2 



2c- 



5(-2 Ui 



c + 



+ 





+ 


+ 


[(■« 


\( 


5 


LVc-i 



1 

:! + C 2 + Cy^ 2 " 
1 3 



12-2c + 



-68j_2_R 

3 + C 



C 4 + 2C 



^+(3i?-2 C -10+^-i 



A 1 A„ „ c-4 
- ! - I 16 + 4c+ — 



28-6i?\ 1 
e 



c-i 



5 + 


"2 


c + 8\ 1" 




2 + C/ c. 




+ 3 4 + C 
2 + C 


n 




+ ( 2C 2 e 2 



Ul 

2 + R 



u 2 



(67) 



where c ~ 0.2274 is another constant found by numerical integration. It may be readily checked that expressions (66) 
and (67) differ from (63) and (64) only by a UV-finite renormalization of the parameters v, u\ and 

In the NP scheme, in contrast with the MS scheme, the renormalized Green functions have an analytic dependence 
on the set of parameters e and A, i.e. they do not have factors of the type ae + bA in denominators. This is in accord 
with the general ideas of the theory of analytic renormalization [11]. 

In the constants Z of the MS scheme with a fixed value of £ = A/e = const the dependence on e is present only in 
the form of poles 1/e, 1/e 2 etc. Contrary to this, in the constants Z of the NP scheme regular terms ~ 1, e, e 2 etc. 
are added to the poles in e. For calculation of the RG functions and the correction exponents u> in Sec. V on rays 
(26) with £ = A/e = const to the order e 2 only terms of order 1/e and 1 are required in the one- loop contributions 
~ 11 to whereas in the two-loop contributions <~ u 2 only terms of order 1/e 2 are 1/e are needed. Expressions (66) 
and (67) are quoted just with this accuracy. 

V. RENORMALIZATION-GROUP REPRESENTATION 

The use of renormalized parameters as such does not solve the main problem of large expansion parameter growing 
with the Reynolds number. It is, however, a necessary step towards the use of the method of the renormalization 
group which allows to solve the problem by effective resummation of the perturbation theory. We shall consider as 
an example the equal-time pair correlation function 



(<Pi{t,*)<Pj(t,x')) = Gij(r) , 



(68) 



14 



which is the most interesting quantity for us in the following. The Fourier transform of this function may be written 
as 

G«(p) = f^(p)G(p), (69) 

where -Py (p) is the transverse projection operator and p = |p|. Dimensional arguments lead to the following repre- 
sentation of the scalar function G(p) from Eq. (69): 

G(p) = v 2 p- d+2 R(s, gi ,g 2 ), 8 =£, (70) 

where R is a dimensionless function of dimcnsionlcss arguments. We want to calculate G(p) in the inertial range of 
the wave number p. Since in the present model (3) the external scale of turbulence has been put equal to infinity, 
this corresponds to the region s = [ijp 3> 1. The perturbation expansion of G(p) contains powers of the parameter s 
whose exponents grow without limit, due to which it is ill-suited for finding the sought asymptotic behavior s — > oo. 
We shall briefly remind solution of this problem within the method of RG. 

Since the fields <& = {ip, ip'} in the present problem are not renormalizcd, the renormalized functions W R differ 
from the unrenormalized ones W = ($...$) only by the choice of variables and the form of perturbation expansion 
(<?i and g 2 instead of gio and 520)1 an d we may write: 

W R (g 1 ,g 2 , v,fj,,...) = W(g 10 , g 20 , v Q , . . . ) . 

Here, e = {fo, gio, 320} is the set of all bare parameters, whereas e = {v, gi, g 2 } are their renormalized analogs, and the 
ellipsis stands for the arguments not affected by renormalization like the coordinates, times etc. The unrenormalized 
functions W do not depend on /1, while the renormalized functions W R do because of the introduction of ji in 
renormalization relations (47). The independence of [i of the functions W is expressed by the equation V^W = 0. 
Here, and henceforth, 2? M = [id^ with fixed bare parameters eo- The equation V^W = written in terms of the 
renormalized functions W R — W and their arguments e, \l is the basic RG equation 

V„W R (g, v,fj,,...) = V RG W R (g, v, . . . ) = 0, (71) 

where T>rc stands for the operation expressed in terms of the renormalized variables: 

V RG = V^ + fod gi + f3 2 d g2 - lv V v , (72) 

where V x = xd x for any variable x. The RG coefficient functions (the anomalous dimensions 7 and the (3 functions) 
in Eq. (72) are defined as 

j a =T>^lnZ a a = {u, gi , g 2 , D 2 } 

t = l,2. (73) 

The term with T> v in Eq. (71) is written with the account of renormalization relation (47) for v and definition 7„ 
(73). From Eq. (73) and renormalization relations (47) it follows 

ft {91 , 92) = gi [-2e - 7 S1 (5i ,92)], (74) 
f3 2 (gi,g 2 )=g 2 [2A--f g2 ( gi ,g 2 )], (75) 
7si = -37^, 7 S2 = 7d 2 - 37„. 

We are interested in the IR asymptotics of small wave vectors p and frequencies u> of the renormalized functions W R 
or, equivalently, large relative distances and time differences in the (t, x) representation [in static objects like (68) - 
(70) dependence on t or u is absent]. It is determined by the IR-stable fixed point g*, at which (3{g*) = for all [3 
functions. The fixed point g* is IR stable, if real parts of all eigenvalues of the matrix = 9/3,/ '^9j\ g - g are strictly 
positive (see, e.g. Refs. [12, 13]). Below it will be shown that in our model (45) the system of two (3 functions (74) 
and (75) in the region of our interest e > 0, A > has an IR-stable fixed point 5* = {gi*,g 2 *} with gi* ^ 0, g 2if ^ 0. 

In its presence it follows from the RG equations (71) that (see, e.g.. Refs. [13, 14]) the sought asymptotics W R 
of the Green function W R has the following property of "IR scaling" [in the (t, x) representation] 



IR 



W fl | Jfl (A- A «t,A- 1 x) - \ A ™W R \ IR (t lX ), 

A W = ^A $ , (76) 
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where x is the set of all coordinate variables and t all times, whereas A > is an arbitrary stretching parameter. 
Summation in expression (76) for A\y goes over all fields $ = {ip, cp'} entering the function W R . In Eq. (76) only 
those arguments of the function W R are explicitly shown which are stretched under given scale transformation. 

The quantities A w and A$ in Eq. (76) are critical dimensions of the frequency to and the fields $ = {ip, ip'}. They 
are all unambiguously (see, e.g.. Refs. [13, 14]) expressed through the quantity 7* = 7„(<7*) - the value of the RG 
function 7^(3), defined in Eq. (73), at the fixed point: 

\ = 1 - it, A v ,=d-A v , 

A ul = 2- 1 * l/ , 7* = 7 „(<fc). (77) 

At the fixed point with gu ^ and 32* 7^ the values 7* = j a {g*) of RG functions (73) are readily found from the 
definition the fixed point 0i(g*) = j3 2 {g*) = together with relations (74) and (75): 7 * x = -2s, 7* 2 = 2A, 7* = 
2s/3, -f* D2 = 2A + 2e. Substitution of 7* = 2e/3 in Eq. (77) leads to formulas (20) and their corollaries (21) for s = 2. 
Thus, in two-charge model (45) with the local renormalization [5] the critical dimensions of the velocity field ip and 
frequency u> at the real value e = 2 retain their Kolmogorov values contrary to the conjecture of the author of Rcf. 
[4]- 

Consider again function (70). It is a particular case of the function W R and satisfies RG equation (71): T>rqG = 0. 
A representation of the solution of Eq. (71) for G(p) convenient for the asymptotic analysis at p — > may be 
obtained with the aid of invariant variables e = e(s, e) corresponding to the complete set of renormalized parameters 
e = {v, (71,32}- They are defined as solutions of the RG equations T>rg£ = with the operator T>rg from Eq. (72) 
and the normalization conditions e = e at s = 1. In terms of the invariant variables the solution of the RG equation 
(71) for G(jp) may be represented as 

G(p) = u 2 p 2 - d R{s, gi ,g 2 ) = D 2 p 2 - d R{l,- gi ,g 2 ). (78) 

The right-hand side of (78) depends on s through the invariant variables e(s, e) only, whose asymptotic behavior in 
the limit s — > 00 - determined by the IR stable fixed point (see below) - is simple: the invariant charges g~\ and g 2 
tend to fixed values g^ = 0(e) and g 2 * = 0(e), whereas the invariant viscosity has simple powerlike asymptotics. It 
may be conveniently determined by expressing the invariant variables e = (v,g~i,g 2 ) in terms of the bare variables 
eo = (^0, ffio, 320) and the wave number p. According to definition, the bare variables eo as well as the invariant 
variables e satisfy the equation Vrq^o = £V e o = 0- The connection between the two sets of parameters is determined 
by the relations 

5io = gip 2e Z gi {g), 

(79) 

valid because both sides in each of them satisfy the RG equation, and because they at s = fj,/p = 1 coincide 
with relations (47) owing to the normalization conditions. Using the connection between renormalization constants 
ZgZl = 1 indicated in Eq. (47) and excluding these constants from the first two relations in (79) we find = 
D10 = giP 2e i^ 3 , and from here 

v=(D 10 p- 2 y-g 1 ) 1 / 3 , 
which for the sought asymptotics s — > 00 with the account of g~\ — > 31* yields 

= (£i /gi*) 1/: V 2e/3 , s^oo. (80) 
Substituting this result in Eq. (78) we obtain 

G{p) ~ (D w /gu) 2/ V- d - 4s/3 R(l, 9*), s - 00. (81) 

This relation will be used in Sec. VI. 

Let us make a remark about relations (79). According to renormalization relations (47), condition D 2 q ~ 320 = 
[see the text following Eq. (44)] imposes the constraint Zjj 2 (g) ~ Z g2 (g) = on the renormalized charges g = {gi, g 2 }. 
From the last relation in Eq. (79) it follows that the invariant charges g = g(s, g) for any value of the variable s = fi/p 
lie on the same constraining surface Z„ 2 = as the initial data g\ s =i — 9- Therefore, the limit values 3* = lim g{s,g) 

lie on the same surface Z g . 2 = 0, i.e. the condition D 20 ~ g 20 = is compatible with the RG analysis. 

All said above is valid for any subtraction scheme, only the explicit form of the RG functions "/ a in Eqs. (73) and 
(74) depends on the choice of the scheme. We shall first quote results of the two-loop calculation in the MS scheme 



u = vZ v (g), 
.920 = 92P~ 2A Z g2 (g) , 
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(Sec. IV) and then briefly discuss modification of formulas in the NP scheme. As said before, no physically significant 
results depend on the choice of the scheme. 

In the MS and MS schemes all RG functions j a are independent of e. In model (45) they depend only on charges 
and the parameter £ = A/e. The two-loop expressions for the constants Z a in Eq. (73) are given by Eqs. (63) and 
(64). In calculation of the quantities 7 a = 2? M lnZ a from Eq. (73) the operation may be replaced by T>rg from 
Eq. (71) and the contributions with and V v omitted, since the quantities Z a do not depend on \i and v. Such a 
calculation yields 

7 „ = 2 (t*i + u 2 ) + 2(4 ^ ± f Ul + 2(5C + 3)ui«2 - 4i? («i + u 2 f + ... , (82) 



_ 2(m+u 2 ) 2 C(13 + 19Qt 
lD2 u 2 (2 + Qu 2 



2(34C+19 + 6C 2 K a 2 



+ — /■ ~ 6u i + ( 13 + 31 0«i«2 

2 + C 



+ 3 

U 2 

Let us remind that u\ ~ <?i and «2 ~ 32 are charges with a more convenient normalization (53), while the ellipsis 
stands for corrections of order 0(u 3 ). 

Substituting quantities (82) and (83) in Eq. (74) we obtain expressions for the (3 functions in the two-loop ap- 
proximation. Then from the conditions (3i(g*) — P 2 (g*) = coordinates of the fixed points g* ~ u* may be found. 
In the framework of the e expansion there are three fixed points [5]: 1) the trivial fixed point m* = 0, u 2 * = 0; 2) 
the "kinetic" fixed point m* = 0, u 2 * ^ 0; and 3) the "Kolmogorov" fixed point ui* 7^ 0, u 2 * ^ 0. In the region 
e > 0, A > of interest for us only the Kolmogorov fixed point is IR stable, for which in the one-loop approximation 

ui»+u 2 * = | + 0(e 2 ), U2* = 9(1 ^ +0(e 2 ). (84) 

From relations (82) and (83) two-loop contributions ~ e 2 to Eq. (84) may be found. We do not quote them, 
because coordinates of a fixed point u* ~ g* do not have direct physical meaning and do depend on the choice of the 
subtraction scheme. Objective quantities independent of the subtraction scheme are the eigenvalues of the matrix 
u>ij — df3i/dgj\ g _ g . In our problem the u) matrix is a 2 x 2 matrix, whose two eigenvalues uj± in the two- loop 
approximation at the Kolmogorov fixed point are 

/ 4 J% 2 - 12C -8\ 2f [4(1 + 30-R - 6 - 12C-9C 2 ] 1 , , , 

w± = C + 7T ± — 7T- ~ e+ - \-3-2R-3C± ^ ]' } e 2 . 85 

V 3 3 / 9 I ^9C 2 - 12C - 8 J V ' 

We quote also for reference the relatively simple expressions for the trace and determinant of the uo matrix, through 
which the eigenvalues u>± are unambiguously expressed: 

Trw = co + +co_ = ^(3C + 4)e- ^(3C + 3 + 2i?)e 2 , (86) 

det u) = uo + u- = ^ (3C + 2)e 2 - ^ (2R + 1)(3( + 2)e 3 . (87) 
o y 

The one-loop contributions ~ e in Eqs. (84) - (86) and ~ e 2 in Eq. (87) were obtained earlier in Ref. [5]. In the 
one-loop approximation this fixed point 5* is IR stable in the sector e > 0, ( > —2/3 in the (e, A) plane. When e > 
and C < —2/3 both eigenvalues (85) are real and have different signs [most easily this may be seen from the one- loop 
contribution in Eqs. (86) and (87)]. With growth of ( upon intersection of the borderline Co = —2/3 both eigenvalues 
become positive and then, upon reaching the next borderline 2(1 — \/3)/3 ~ —0.488, the argument of the root in Eq. 
(85) becomes negative, i.e. the fixed point becomes an IR-attractive focus with u>± — a±ib with a > 0. It remains 
such until the next borderline 2(1 + v / 3)/3 ~ 1.821 is reached, upon passing which the root argument in Eq. (85) 
becomes positive again and both eigenvalues u>± real and positive. For our "physical" ray £ = 1/4 (d = 3) the fixed 
point is an IR-attractive focus. 

What was said above refers to the one-loop approximation. The account of the two-loop corrections in Eqs. (85) 
- (87) leads to a deformation of the borderlines of the region of IR stability, but the "physical" segment of ray (26) 
with C = 1/4, 0<e<2 still remains in this region. 
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Consider now renormalization in the NP scheme, in which the constants Z are determined from the normalization 
condition (65). The two- loop expressions for the constants Z with the necessary accuracy [see the text following Eq. 
(67)] are given in Eqs. (66) and (67). With the use of them together with definitions (73) we obtain the following 
expressions for j a : 

7„ = (2 + 3A-ce)ui + (2 + cA +3A)n 2 -4 ( Ul + u 2 f {2R + 1) , (88) 



7r> 2 



2 + (7 - c) A + (4 - 2 c) e 



u 2 



+ 2 



2 + 5A + (2 - c)e 



Ml + 



2 + (3 + c)A 



«2 



_ 4 ( Ul W(2fl+l) 

U 2 



where the notation is the same as in Eqs. 66) and (67). The RG functions (88) and (89), contrary to their analogs in 
the MS scheme, do not contain factors like ( + const in denominators, i.e. they are analytic in the pair of parameters 
e, A, which is a consequence of similar analyticity of the renormalized Green functions. Coordinates of fixed points 
~ g* obtained from Eqs. (88) and (89) in the one- loop approximation keep the form of Eq. (84), but the two- loop 
contributions (which we do not quote) differ from analogous contributions in the MS scheme. The eigenvalues to± of 
the matrix u>, however, remain exactly the same as in the MS scheme, because these quantities do not depend on the 
subtraction scheme. 

In conclusion, we note that in an attempt to use the NP scheme (43) in the model [4] with nonlocal renormalization, 
the inconsistency of this model in terms of the RG functions -f a would show in the form of poles 1/e in the two- loop 
contributions. 



VI. SKEWNESS FACTOR AND KOLMOGOROV CONSTANT 



The exponent of the power of the wave number in Eq. (81) is determined exactly and does not have corrections 
in the form of higher powers of e. At the physical value e = 2 this exponent assumes the Kolmogorov value. To 
find the Kolmogorov constant, the amplitude of this function has to be calculated, which, however, can be done only 
approximately, because the corresponding e expansion does not terminate. In calculation of the amplitude, apart 
from technical difficulties at two-loop order, a principal problem arises as well. It is connected with the necessity 
to express the answer for G(k) in terms of the energy injection rate £ instead of the parameter D\q of the forcing 
correlation (44). The connection between D\$ and £ is determined by an exact relation expressing £ in terms of the 
function d/(k) in the correlation function (2) 

£=S4r fdkdf{k). (90) 



2(27r) fl 

Substituting here function (44) with Z?2o = [see the text following Eq. (44)] and introducing the UV cutoff 
k < A = (5/Vq) 1 / 4 (the inverse dissipation length), we obtain the following connection between the parameters £ and 
D w 

«p_-«)A-« 

Idealized injection by infinitely large eddies corresponds to df(k) cx S(k). More precisely, according to Eq. (90) 



2(2Tr) d £S(k) 

«/W = 

In view of the relation 



d/(fc)= V '"j_7 W - (92) 



S(k) = lim (2n)- d [ d*(Ax) 2e - 4 exp(ikx) = S^k~ d lim [(4 - 2e)(fc/A) 4 - 2e ] , 

e^2 J ' e— »2 

the powerlike injection with df — Dmk A ~ d ~ 2e and the amplitude Dio from Eq. (91) in the limit e — > 2 from the the 
region < e < 2 gives rise to the 5 sequence (92). 
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Relation (91) reveals that at fixed £ the quantity D\$ depends on e and it is necessary to take this dependence into 
account in the construction of the e expansion for the Kolmogorov constant. On the other hand it shows that the 
quantity R(l,g*) from (81) must have a singularity of the form (2 — e)~ 2 / 3 in the limit e — > 2: only in this case the 
Kolmogorov constant in the model with the injection df — D w k 4 ~ d ~ 2e and the amplitude D w from Eq. (91) shall 
have a finite value in the limit e — > 2. The measurable experimental Kolmogorov constant Ck in terms of the model 
with such pumping corresponds to the limiting value e = 2, and we want to define its generalization Ck(e) for the 
whole interval < e < 2. Obviously, such a generalization cannot be done unambiguously, because it is not possible 
to define unambiguous dependence of the parameter D w in Eq. (91) on e at a fixed value of £. 

Let us explain this in more detail. When deriving relation (91) we assumed that integral (90) for the injection 
df = D w k 4 ~ d ~ 2e has an upper cutoff equal to the inverse dissipative length A = (£ /i/q) 1 / 4 . Such a cutoff is natural, 
but at the same time only orders of magnitude may be discussed, of course, not the exact values. Therefore, there 
is nothing to prevent to replace in Eq. (91) the cutoff parameter A by a A with a coefficient a of the order of unity, 
which yields the extra factor a 2e ~ 4 on the right-hand side of Eq. (91). This factor tends to unity at e — > 2, hence 
it does not affect the the physical (real) value of the Kolmogorov constant Ck(s = 2), but it does affect coefficients 
of the hypothetical e expansion of the function Cr-(e). Generalizing these observations it may be stated that the 
physical content of the theory is not changed, if to the right-hand side of Eq. (91) an extra factor F(e) with F(2) = 1 
is added. In Ref. [18] (see also [13, 14, 17]) relation (91) without the extra factor F(e) was regarded as the definition 
of the quantity D w . Other approaches to the definition of the function Ck(s) and its e expansion [19] - [25] may be 
reduced to the introduction of a particular function F(e) with F{2) = 1 on the right-hand side of relation (91). 

Thus, e expansion of the Kolmogorov constant in the model with the powerlike injection is not defined unambigu- 
ously. However, physical quantities independent of the amplitude D w (universal quantities) do have a well-defined e 
expansion. The skewness factor 

S = 5 3 /5 2 3/2 , (93) 
is an example of such a quantity. In Eq. (93) S n are structure functions defined by relations 

S n (r) = ([p r (t,x + r)- Vr (t,x)] n ), Vr =^^.. (94) 

According to Kolmogorov theory, the structure function S^r) in the inertial range is of the form 

S 2 (r) = C K t'\ 2 l\ (95) 

where Ck is the Kolmogorov constant with a simple connection with the Kolmogorov constant of the energy spectrum 
[16]. Since the experimental evidence for anomalous scaling [i.e. deviation of the power of r from the Kolmogorov 
value | in Eq. (95) in the inertial range] in the structure function S 2 (r) is still controversial and in any case this 
deviation is small [26], we shall use the Kolmogorov asymptotic expression (95) in the following analysis. 
The structure function Ss(r) may be found exactly in the inertial range [16]: 

s >u = -wh) Ir > (96) 

which allows - with the account of (93) and (95) - to relate the Kolmogorov constant and the skewness factor 
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C K = 



d(d + 2)S 



2/3 

(97) 



Of the three quantities ^(r), Ss(r) and S only S has a unique e expansion. Thus, relation (97) (valid only for the 
physical values = 2) might be used to determine Ck using the calculated value S(e = 2). 

To find the RG representation of skewness factor (93) it is necessary to have RG representations of the functions 
52 (r) and Ss(r). The function 5*2 (r) is connected with the Fourier transform of the pair correlation function G(k) by 
the relation 



/<7k 



1 (k-r) 



(At) 



1 {l-exp[i(k •!■)]}, (98) 



therefore its RG representation may be found on the basis of RG representation (81). An analogous RG representation 
in the inertial interval may be written for the function Ss(r). It is more convenient, however, to use the following 
exact result, an analog of expression (96): 



3(d-l)r(2-e) (r/2) 2s - 3 D 
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This relation is a manifest demonstration that the amplitude of the structure function, expressed in terms of D\q, 
has a singularity at e — > 2, in this case it is <~ (2 — e) _1 . On substitution of Eq. (91) in Eq. (99) this singularity 
cancels the corresponding zero on the right-hand side of Eq. (91), leading for Ss(r) to an expression finite at e = 2 
and coinciding with Eq. (96). 

Relations (81), (98) and (99) might serve as the basis for construction of the e expansion of the skewness factor 
(93). However, an additional difficulty arises on this way. The point is that the powerlike dependence £2 (r) ~ ?- 2_2£ / 3 , 
determined from Eqs. (81) and (98), is only valid when e > 3/2, because for e < 3/2 integral (98) diverges at k — > 00 
[this means that the main contribution to 52 (r) in this case is given by the term (<£ 2 (t, x)) independent of r}. However, 
the derivative rd r S2(r) is free from this flaw, because, according to Eq. (98), 



rd r S 2 (r) = 2 



(27T) 



(k • r) 2 



{krf 



(k • r) sin(k • r). 



(100) 



Integral (100) is convergent for all < s < 2. On the other hand, at the physical value e — 2 the amplitudes in 52 (r) 
and rd r S2(r) differ by a trivial factor |, therefore in Refs. [1-3] for the construction of the e expansion the following 
analog of the skewness factor was used 



Q(e) = r9rS2{r) rdrS2(r) 



|5 3 (r)| 2 /3 (-5 3 (r)) 2 / 3 ' 



(101) 



The Kolmogorov constant and the skewness factor are expressed through the value Q(e = 2) according to Eqs. (93), 
(95) and (96) by the relations 



3Q(2) 
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d(d + 2) 



2/3 



3Q(2) 



3/2 



(102) 



Quantity (101) may be calculated both in the double (e, A) expansion and in the usual e expansion. In the former 
case the corresponding expansion is obtained on the basis of relations (81), (99) and (100) in the form: 



Q( £ ,C)=e 1/3 5> fe(C)£ *. 

fc=0 

The usual e expansion of the quantity Q for dimensions d > 2 has been obtained in Ref. [1]: 

00 

Q(e,d)=e 1 ^Q k (d)e k . 



(103) 



(104) 



fc=0 



The connection between expansions (103) and (104) is revealed by investigation of singularities of the coefficients 
Qk(d) in Eq. (104) at d — > 2. An analysis of these singularities shows that in the vicinity of d — 2 = 2 A = these 
coefficients may be expressed in a Laurent expansion 

00 

Q fc (d) =]r ft; A'- fe . (105) 

(=0 

Substitution of expression (105) in Eq. (104) leads to the representation 

00 00 

Q( £ ,rf)= £ 1 /3^^ (e/A) fe toA i. (106) 

k=0 1=0 

Changing variables in Eq. (106) to e and ( — A/e, we arrive at expansion (103), in which 

00 

1=0 



Relations (105) and (107) show that the alternative e expansions (103) and (104) sum different infinite subsequences 
of double sum (106). In Ref. [3] a procedure of improvement of the e expansion was proposed with the use of the 
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1 k 

FIG. 2: Summations in the calculation of Qj}f in Eq. ( 109). Terms in the double sum (106) taken into account in Q ™^ and 

Qi , correspond to the dashed horizontal and and vertical stripes, respectively. The correction term <5Q' n ' corresponds to sum 
over the double-dashed square. 

mutually complementary information about the quantity Q contained in the partial sums of expansions (103) and 
(104) 

7i — l n—1 

q ( :X = £l/3 E **(0e*> Qi n) = s 1/3 E w)^ fe , ( 108 ) 

fc=0 fe=0 

where n > 1 is the number of loops. 

Terms in the double sum (106) taken into account in ano ^ have been schematically plotted in Fig. 2 in 

the form of dashed horizontal and and vertical stripes, respectively. 

All terms in the dashed area will be taken into account in the effective quantity 

Q^ f = Q^+Q^l-SQ^\ (109) 

where 

n—1 n—1 

aQ(»W/3££ (e/A) fc gfclA « 

fc=0 1=0 

is a subtraction term necessary to avoid double counting of terms with k < n — 1, I < n— 1 (the double-dashed square 
in Fig. 2). It may be found by taking the corresponding number of terms from expansions (105) or (107). From the 
point of view of the usual e expansion (104) relation (109) may be interpreted as follows: in the n—1 first terms 
of the expansion the coefficients Qk{d) from Eq. (104) are calculated exactly, but in all higher-order terms (k > n) 
approximately with the account of n — 1 first terms of their Laurent expansion (105). 

Our two-loop calculation of the e, A expansion of the quantity Q together with the two-loop calculation of Ref. 
[1] allowed to obtain an improved e expansion of the quantity Q at second order of perturbation theory [3]. For the 
Kolmogorov constant calculated according to Eq. (102) for d = 3 it led to the result quoted in Table I. 

In Table I we have quoted for comparison the values of the Kolmogorov constant calculated according to Eq. (102) 
at first and second order of the usual e expansion (C £ ), the double e, A expansion (C e .a), the contribution Cs in 
Eq. (102) from the correction SQ n in Eq. (109) and the value C e // obtained from relations (102) and (109). In all 
the cases quoted the recommended experimental value of the Kolmogorov constant C exp = 2.01 [27] lies between the 
values of the first and second approximation. However, the difference between these values is rather significant both 



2-loop 
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TABLE I: One and two-loop values of the Kolmogorov constant in the usual e expansion (C e ) and the double e, A expansion 
(C £ ,a); the contribution Cg in Eq. (102) from the correction 8Q (n> in Eq. (109), and the value C eff from Eqs. (102), (109). 



n 


c € 


C e ,A 


C s 




1 


1.47 


1.68 


1.37 


1.79 


2 


3.02 


3.57 


4.22 


2.37 



in the e expansion and in the (s, A) expansion, let alone the leading terms of the s expansion of the latter. For the 
improved e expansion, i.e. for the quantity C e ff = C e + C £: a — Cs calculated according to Eqs. (109) and (102), 
however, this difference is about three times smaller leading to far better agreement with the experimental data. 

VII. CONCLUSION 

In conclusion, we have presented a detailed comparison of two different renormalization schemes for the stochastic 
Navier-Stokes problem near two dimensions. By explicit two-loop calculation we have shown that the nonlocal scheme 
of Ref. [4] cannot consistently be carried out beyond the leading one-loop approximation. On the contrary, our two- 
loop results confirm the consistency of the local renormalization scheme of Ref. [5] based on the general principles of 
the theory of UV renormalization. 

The detailed explicit two-loop analysis of different renormalization schemes presented here is all the more important, 
because the inconsistent renormalization of nonlocal terms in dynamic models continues to appear in the literature 
[28,29]. 

The correct choice of the renormalization scheme is vital for a proper account of the effect on structure functions 
of the additional singularities appearing in the field-theoretic model in the limit d —> 2. Using the consistent local 
renormalization scheme, we have shown that a proper account of the " nearest singularity" in the coefficients of the s 
expansion (104) leads to a significant improvement of the results of the two-loop RG calculation at d = 3. We have 
analysed the effect ot this procedure at other d as well. It turned out to reduce significantly the relative contribution 
of the two- loop correction in the whole range considered oo > d > 2.5. At the same time this contribution remained 
large at d = 2, which we think to be an effect of singularities at the next exceptional dimension d = 1. 

The proposed procedure of approximate summation of the e expansion is, of course, applicable not only to the 
calculation of Q(s), but all universal quantities such as dimensions of composite operators. 
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APPENDIX: A RENORMALIZATION AND (e, A) EXPANSION ABOVE TWO DIMENSIONS 

As it was explained in Sec. IV, in the two-charge model (45) in some graphs the wave-vector integrals diverge at 
large wave numbers. To regularize such integrals, it is necessary to introduce a cutoff parameter A. This may be 
done, e.g., by restricting Fourier components of the velocity field <p to wave numbers less than A in functional (45), 
which automatically brings about the corresponding sharp wave-vector cutoff in the bare response function (5) and 
in the bare correlation function (46). It was already explained in Sec. IV that all such A divergences are "nearly 
logarithmic" and appear in the results in the form of powers A" with small (of the order of e for A/e = const) positive 
exponents a. 

The elimination of the A divergences may be reduced to a renormalization of the bare parameters. Denoting for 
brevity the whole set of parameters by e we introduce the notion of "primary bare parameters" eo = {j>o,£>io = 
ffio^o' * = 1 5 2} and "secondary bare parameters" eo = {vq,Diq = 9i0^q, i = 1,2} (see Ref. [13]). The original model 
is defined by a functional of the type of Eq. (45) with the A cutoff introduced and with the "primary bare parameters" 
e : 

= <p'(D w k 2 - 2A - 2s + A>ofc V/2 + <p'[-dt<p + i>od 2 <fi - (<pd)<p] . (A.l) 
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Renormalization of this model may be carried out in two steps: the first is the A renormalization with the aim of 
removal of all A divergences. This amounts to a reorganization of the bare parameters e — > e , in which the secondary 
set of parameters is expressed as functions of the primary set: eo = eo(eo,A); and vice versa eo = eo(eo,A). The 
correspondence between the two sets eo and eo is bijective perturbatively, therefore any of them may be chosen as the 
set of independent variables. 

If the parameters eo are chosen as independent, then in the Green functions V of model (A.l) expressed in terms of 
e and A there will be no A divergences left (they all will be concentrated in the formulas connecting eo and eo) and 
the limit A — > oo may be taken in them with the result of eliminating the cutoff parameter A completely from the 
theory Trace of the UV divergences which brought about the positive powers of A remains, however, in the form of 
singularities in e in the A-renormalized quantities. This happens because in the A renormalization only terms strictly 
growing as powers of A are removed and collected in the renormalization constants. These terms contain singularities 
in e and A, although the unrenormalizcd quantities with fixed A were regular functions of e and A. Consequently, in 
the A-renormalized quantities there must be terms left which are singular in e and A, but remain finite in the limit 
A — > oo. Thus, the A renormalization is a way to trade UV divergences in the form of positive powers of the UV 
cutoff A for poles in e, A and their linear combinations in such a way that in the A-renormalized quantities the limit 
A — > oo may be taken. 

The basic conjecture is that the results obtained in this manner for the graphs of the Green functions r(eo, A = oo, ...) 
(the ellipsis stands for the rest of the arguments like frequencies and wave vectors) are exactly the same as those 
obtained in the "formal scheme" , i.e. by analytic continuation of all integrals without A divergences on the parameter 
A from the region of small A < (more accurately — 2e < A < 0). In this scheme, the unrenormalizcd action is 
functional (45). Such an analytic continuation might be carried out without any reference to the model regularized with 
the explicit wave-number cutoff A, which is common practice in field theories of particle physics. There, however, 
it is the renormalized parameters which are the physical ones and their bare counterparts together with the UV 
cutoff unphysical auxiliary quantities. In our case unrenormalized parameters are the physical ones and therefore 
it is important, in principle, to keep track on their relation to the (auxiliary) renormalized parameters, because the 
fixed-point values of the latter remain in the asymptotic expressions for various correlation functions and the like. 

The next step after the A renormalization is the e renormalization with the goal of removal from all Green functions 
T(eo, A = oo, ...) poles in e for A/e = const. It is carried out by the transition from the "secondary bare parameters" 
eo (the same notation was used in Sec. IV) to the renormalized parameters e = {v, gi, 52} according to relations (47). 

The procedure of the e renormalization was discussed thoroughly in Sec. IV. Let us now explain in more detail 
the procedure of the A renormalization: the transition from the primary bare parameters eo to the secondary bare 
parameters eo- We emphasize that at this stage we are interested in the A divergences only and regard e and A as 
fixed parameters without any investigation of singularities in these parameters tending to zero. We shall consider 
the parameters eo in the graphs of the functions T of model (A.l) expressed in terms of eo and A through the 
renormalization relations 



(A.2) 





= 510^0 = 


510^0 — ^10 ! 




= .920 = 


92oVqZd 2 — D20ZD2 7 


510 


= 9ioZ gi , 


ff20 = ff20 Zg 2 , 


^0 


= VqZ v , 


ZgiZ v — 1, Z g2 Zl = 



^D 2 



similar to Eq. (47). The dimensionless renormalization constants Z in Eq. (A.2) are functions of e and A expressed 
in the form of series in D i0 <~ g i0 . The corresponding dimensionless expansion parameters are the following analogs 
ofEq. (48): 

DioSd _ D20S11 , . „s 

ai =32^' a2 = 32*gA-*A • (A ' 3) 

Therefore, the constants Z in Eq. (A.2) assume the form 

z„,d 2 = i + J2 cSr^r^ 2 (a.4) 

711 >0, 712>1 

with the dimensionless coefficients ^ depending on e and A/e = £ only (but in a singular manner!). In 
expansion (A.4) not all possible terms are included, but only those which are "A divergent", i.e. those with a positive 
power of A in the product a™ 1 a!? 2 . From Eq. (A. 3) it follows 

d$ 1 a%*~A- a , a = 2(n 2 A-n 1 e), (A.5) 
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therefore, for e > 0, A > 0, in the A-divergent terms with a > in Eq. (A. 4) the inequality ri2 > 1 holds, i.e. at least 
one factor with D 2 o ~ D 2 o from (A.l) is present. 

From this it follows, in particular, that to the real value -D20 — in Eq. (A.l) [see the text after Eq. (44)] it 
corresponds Z? 2 o = in Eq. (45), which justifies the derivation of Eq. (91) from Eq. (90) in model (45). We also note 
that the operation in Eq. (73), defined in Sec. V as 2? M = fid^ with fixed parameters e , in terms of model (A.l) 
has to be understood as [id^ with fixed eo and A. These definitions arc equivalent, because the parameter fj, does not 
enter in renormalization relations (A. 2). 

For the A renormalization (A. 2) analogs of relations (49)-(52) may be written and the corresponding Z calculated 
at two-loop order. We shall not quote the corresponding results, because explicit expressions connecting the primary 
(eo) and secondary (eo) bare parameters is unimportant for the RG analysis of the IR asymptotic behavior in Sec. 
V, which is carried out in terms of bare parameters eo and renormalized parameters e. 

In the NP scheme the normalization condition (65) may be imposed in the A-renormalized model in the same way 
as just described for the MS (or MS) scheme. It is not difficult to see, however, that in the NP scheme the very 
procedure of the A renormalization is actually not necessary. The point is that in this scheme from the quantity to be 
renormalized its value at the normalization point is subtracted which automatically leads to a quantity without any 
UV divergences and thus with a finite - and regular in e and A - limit, when A — > 00. For renormalized correlation 
functions the result is the same as after A renormalization, subsequent limit A — > 00 and final renormalization in 
the NP scheme. Therefore, the RG functions 7 and j3 are also the same, since their expressions in terms of the 
renormalized correlation functions coincide in both cases. 
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